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ABSTRACT 

This  thesis  describes  an  alternating  renewal  process  and 
a  computer  program  that  simulates  a  transient  alternating 
renewal  process.   The  purpose  of  the  simulation  is  to  test 
analytical  approximations  for  the  efficiency  of  a  system 
described  by  such  a  process. 
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I.   DESCRIPTION  OF  AN  ALTERNATING  RENEWAL  PROCESS 

This  thesis  describes  an  alternating  renewal  process  and 
gives  a  computer  program  that  simulates  a  transient  alter- 
nating renewal  process.   Consider  a  system  that  can  be  in 
one  of  two  states,  up  or  down.   The  system  starts  up  and 
alternates  between  up  and  down,  remaining  in  each  state  a 
random  time.   The  successive  times  in  each  state  are  jointly 
independent.   The  sequence  of  up  times  are  identically 

distributed  with  distribution  F   and  the  sequence  of  down 

u  M 

times  are  identically  distributed  with  distribution 

F,.   X.  and  Y.  are  the  corresponding  i    duration  of  time 

the  system  is  in  the  up  state  and  down  state.   The  means  of 

the  two  distributions  are  assumed  to  be  finite  and  are 

denoted  by  y   and  y  ,  respectively.   The  up  rate  and  down 

rate  are  denoted  by  A   and  A-,,  and  are  related  to  the 

J       u       d 

means  as 


y   =  r=-  =  mean  time  up 

Mu    A  * 

u 

y ,  =  —-   =  mean  time  down 
d 


(1) 
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Figure  1.   An  Alternating  Renewal  Process. 


Two  quantities  of  interest  concerning  an  alternating 

renewal  process  are  the  availability  and  the  efficiency. 

The  availability,  A(t),  is  the  probability  that  the  system 

will  be  up  at  time  t.   The  efficiency,  E(t) ,  is  the  ratio 

of  the  expected  amount  of  up  time  during  [0,t]  to  t.   The 

relationship  between  E(t)  and  A(t),  which  will  be  derived 

later,  is 

t 
/  A(x)  dx 

E(t)  =  ^ (2) 

t 

An  alternating  renewal  process  has  two  imbedded  renewal 

processes.   Reference  1  gives  a  detailed  description  of  a 

renewal  process.   One  process,  call  it  a  type  1  process,  has 

as  its  renewals  the  events  when  the  system  goes  from  the  up 

to  the  down  state.   This  is  a  general  renewal  process  with 

the  first  interevent  distribution  F,  =  F   and  the  remaining 

1     u  * 

interevent  distributions  all  being  F  =  F   o  F , ,  where  'o1 

J       u    d 

denotes  convultion.   The  duration  of  the  i    interevent  time 

is  Y.  n  +  X.  for  i  >  2.   A  second  imbedded  renewal  process 
l-l     l       —  * 

has  as  its  renewals  the  events  when  the  system  goes  up.   The 

common  interevent  distribution  is  F  =  F   o  F,.   This  will 

u    d 

be  denoted  as  a  type  2  process. 

The  expected  number  of  renewals  during  [0,t]  is  denoted 
by  M(t)  and  is  called  the  renewal  function.   If  F  is  a 
function  of  bounded  variation,  that  is,  the  difference  of 
two  monotone  functions,  and  vanishes  to  the  left  of  zero, 
then 
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Figure  2.   Type  1  and  Type  2  Renewals 
and  the  Interevent  Times. 


t 


F  (s)  =  /   exp (-sx)dF(x) 
0- 

is  defined  as  the  Laplace-Stiel jes  Transform  of  F.   If  the 

times  between  renewals  are  distributed  as  F,,F,...,  then 


Fl(s) 
M  (s)  = 


1  -  F*  (s)  * 


(3) 


For  a  derivation  of  this  expression  see  Ref.  1. 

For  the  type  1  imbedded  renewal  process  with  renewal 


function  M, (t) , 


M*(s)  = 


F*(s) 
u 


1  -  F*(s)  F*(s) 
u     d 


For  the  type  2  imbedded  renewal  process  with  renewal 

function  M  (t) , 

F*(s)  F*(s) 
»**  /  \       u     d 
M  (s)  = . 

'  1  -  F  (s)  F*(s) 

u     d 

Define  the  random  variable  S (t)  such  that 

1   if  the  system  is  up  at  t 
S(t)  = 

.  2   if  the  system  is  down  at  t. 
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Figure  3.   Alternating  Renewal  Process  with  Two 

Imbedded  Renewal  Processes. 


Let  N   =  number  of  renewals  of  type  i,  i  =  1,2,  that  occur 
in  [0,t] .   Then 


S  (t)  =  1C  -  NT;  +  1. 


If  the  expectation  is  taken,  then 


E[S(t)]  =  E[n£]  -  E[n£]  +  1 


=  M2(t)  -  Mx(t)  +  1  =  A(t) 


The  Laplace-Stiel jes  Transform  of  (4)  gives 


A*(s)  =  M*(s)  -  M*(s)  +  1 


(4) 


F*(s)  F*(s) 
u      d 

1  -  F*(s)  F*(s) 
u      d 


F*(s) 
u 


1  -  F*(s)  F*(s) 
u     d 


+  1 


A*(s) 


1  -  F  (s) 
u 


1  -  F*(s)  F*(s) 


(5) 


Let  the  random  variable  Q(t)  be  the  amount  of  up  time 
during  [0,t]  and  R(t)  be  the  expected  amount  of  up  time  over 
the  interval  [0,t].   Then 


Q(t)  =  /   S(x)  dx  . 
0 

If  the  expectation  is  taken,  then 


R(t)  =  E[Q(t)  ]  =  E[/   S(x)dx] 

0 


=  /   A(x)  dx 
0 

from  which  (2)  immediately  follows 


/   A(x)  dx 
E(t)  =  5 . =  R(t)/t  . 


If  (5)  can  be  inverted,  then  availability  and,  hence 
efficiency,  can  be  obtained.   For  example,  if  the  distribu- 
tions of  the  up  and  down  times  are  exponential  with  cor- 
responding rates  A   and  A,,  then 


Fu(t)    =    1    -   exp(-Aut) 


Fd(t)    =    1    -   exp(-Adt) 


Fr(s)    -    A-T-i- 


i   =   u,d 


(6) 


gives 


1    - 


A"(s)    = 


u 


A    +s 
u 


\VS)U+S 


1- 


A  ,    +    s 
d 


A      +    A  ,    +    s 
u  d 


Inverting  yields 


A(t)    = 


K^^lh-'v^1] 


Use   of    (2)    gives 


Eft)    = 


lXd+ 


^  +  (vt)(exp[-a^)x,l. 


dx 


u 


Au+Xd        t(X    +AJ 
u      d 


2    (1-  exp[-(Xu+Xd)t]) 


A  (s)  is  not  always  easy  to  invert.   For  instance,  if 
the  down  distribution  is  k-erlang,  then 

/  ^    k 


1  - 


u 


A*(s)  = 


X  +s 
u 


(vs)  (v 


,.,, 


(X  +s)  (x  ,+s) -XX 
\  u  /  *  d  /      u  ( 


(7) 


The  inversion  of  (7)  can  be  obtained  in  principle  but  it  is 
not  easily  expressed  in  closed  form.   If  the  down  distribu- 
tion is  log  normal,  a  reasonable  assumption  for  many 
systems,  then  A(t)  cannot  be  obtained  in  closed  form. 

If  the  system  has  been  operating  for  a  long  period  of 
time  it  will  settle  down  into  a  steady  state.   For  large  t 
an  asymptotic  approximation  for  E(t)  can  be  found  using 


/  A(x)  dx 
lim   0 

t  -*■  °°  +- 


1   lim 


0+ 


s    A  (s) 


and 


d  F* (s) 


d  s 


=  -  y  =  F*«  (s) 


s  =  0 


,  F is  any  distribution 
s=  0 


Letting  n  =  1,  then 


/  A(x)  dx 
lim  0 1_  lim 

t  -*■  °°     +.        1!  s  ■*  0+ 


1  -  F  (s) 

u 


1  -  F*(s)  F*(s)J 
u      d 


Use  of  L1 Hospital's  rule  for  limits  in  the  indeterminate 
form  gives 


/   A(x)  dx 
lim   C) _  lim 

t  -*  oo  '       4-  S  •*■   0  + 


-F  '(s) 
u 


-[F*«  (s)F*(s)  +F*(s)F*'  (s)] 


u 


Ku   d 


(8) 


The  derivation  of  (8)  was  independent  of  the  type  of  distri- 
butions.  If  the  system  is  in  a  transient  state,  and  equi- 
librium has  not  been  reached,  then  the  asymptotic 
approximation  is  not  valid  and  other  approximations  must  be 
made.   Hence,  a  need  exists  for  testing  the  approximations. 
The  program  described  in  the  next  section  simulates  an 
alternating  renewal  process  for  this  purpose. 
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II.   A  SIMULATION  OF  AN  ALTERNATING  RENEWAL  PROCESS 

This  section  describes  a  simulation  of  an  alternating 
renewal  process.   Three  distributions  are  used. 
Exponential 

F(t)  =  1  -  exp(-At)     t  >  0 

=  0  otherwise. 


k-erlang 


k-1  (St)1  exp(-Bt) 


F(t)  =  1  -   I   - t  >  0 


i=0 


i! 


=  0  otherwise 

Log  normal 

F(t)  =  /   \==:     exp[ -1—  (£n  x-a)2]  dx    t  >  0 

0   ax/27T         2a 

=  0  otherwise. 

Two  combinations  of  the  distributions  are  used.   In  one 
combination  the  up  times  are  distributed  exponentially  and 
the  down  times  are  distributed  k-erlang.   In  the  other  sit- 
uation, the  up  distribution  is  the  same  but  the  down 
distribution  is  log  normal. 

R(t) ,  the  expected  amount  of  up  time,  and  E(t) ,  the 
efficiency  of  the  system  over  the  interval  [0,t],  are  the 
output  parameters  and  are  provided  for  ten  time  intervals. 
Analytical  approximations  for  E(t)  are  given  for  comparison. 


11 


The  input  parameters  vary  with  the  choice  of  down 
distribution.   There  are  some  inputs  common  to  both  situa- 
tions:  At  the  unit  time  interval,  X      the  up  time  rate, 

X-,  the  down  time  rate,  the  number  of  simulations  and  runs 
d 

desired.   If  the  k-erlang  distribution  is  selected  the 
parameter  k  must  be  provided.   If  the  log  normal  distribu- 
tion is  used  the  variance  must  be  given.   The  computer 
program  gives  more  information  on  the  inputs. 

The  program  functions  generally  as  outlined  in  the  flow- 
chart.  Each  simulation  has  these  basic  steps.   First  the 
up  cycle  is  entered.   While  in  the  up  cycle  and  upon 
exiting,   checks  are  made  to  determine  if  a  time  period  has 
been  spanned.   If  not,  the  next  step  is  the  down  cycle. 
While  in  the  down  cycle  and  upon  exiting,  tests  are  made 
again  to  determine  whether  or  not  the  end  of  a  time  period 
has  been  reached.   If  not,  the  program  returns  to  the  up 
cycle.   If  at  any  check  the  time  interval  over  which  the 
simulation  is  to  occur  has  been  spanned,  and,  if  the  desired 
number  of  simulations  has  not  been  completed,  then  another 
simulation  begins.   After  the  required  number  of  simulations 
have  been  completed,  the  output  parameters  are  computed  and 
printed. 

The  program  was  written  in  FORTRAN  IV  and  used  on  an 
IBM  360/67  computer. 
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